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Importance  of  Collision  Physics 


Important  Collisions  in  Spacecraft  Propulsion: 


o  Discharge  and  Breakdown  in  FRC 
o  Collisional  Radiative  Cooling/Ionization 
©  Combustion  Chemistry 

Common  Features  in  Spacecraft  Collisions: 

©  Relevant  Densities  Spanning 
Many  Orders  of  Magnitude  —  6+ 

©  Transitions  from  Collisional  to  Collisionless 

©  Tiny  Early  e~  or  Radical  Populations  Critical 
to  Induction  Delay 

©  Many  types  of  Inelastic  Collisions  with 
Unknown  Effects  on  Distribution  Shapes 
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Need  Low  Noise  &  High  Dynamic  Range 
Collision  Algorithms 


Standard  Collision  Models 


Previous  Collision  Methods: 

o  Monte  Carlo  Collisions  (MCC) 

-  Particles  Collide  with  Background  “Fluid” 

-  Often  Used  in  Plasma/PIC  Simulation 

-  lon-e~  Collisions  Assume  Stationary  Ions 

-  No  Conservation/Detailed  Balance 

o  Direct  Simulation  Monte  Carlo  Collisions  (DSMC) 

-  Most  Modem  Versions  use  No-Time  Counter  (NTC)  Method 

-  Conservative/Reversible  Collision 

-  Satisfies  Detailed  Balance 

-  Subset  of  Possible  Collisions  Sampled 

-  Random  Selection  vs  Zy  for  All/Nothing  Collision 


All  Random  Flip  vs  Number  of  Collisions:  Zy  =  ^  (crv)  dt 
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Variable  Weights  for  Dynamic  Range 


Continuum  to  Discrete  Representation: 

o  Many  Particles  ^  Continuous  Distribution 
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Variable  Weights  for  Dynamic  Range 


Continuum  to  Discrete  Representation: 


o  Many  Particles  — )>  Continuous  Distribution 

o  Discretized  VDF  Yields  Vlasov 
But  Collision  Integral  Still  a  Problem 

©  Particle  Methods  VDF  to  Delta  Function  Set 

©  Collisions  between  Discrete  Velocities 

©  But  Poorly  Resolved  Tail 

(Tail  Critical  to  Inelastic  Collisions) 

©  Variable  Weights  Permit  Extra  DOF  in  Tails 

Variable  Weight  “All-or-Nothing”  Collisions? 

Physically  Inconsistent! 

(Mixing  Violates  Momentum/Energy  Conservation) 
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Review  of  Fractional 


NTC  Collisions: 

o  (Collision  Rate  Volume): (Cell  Volume) 
o  Select  Fraction  of  \ N 2  Possible 


Fractional-NTC  Collisions: 


Collisions 


Zij=^(av)ljdt=^(*v)ljdt 

Pij  =  w  (av)ij  dt  IV ceu 
Pmax  =  W  (av)™  dt/Vcell 

N select  =  ^Fniav)™  dt/Vcell 


s  > 


:  □ 


'O  c*s  (y 


Review  of  Fractional  Collisions 


NTC  Collisions: 


o  (Collision  Rate  Volume): (Cell  Volume) 
o  Select  Fraction  of  |V2  Possible 
o  Probability  of  Event  Ratio  of  Volumes 


Fractional-NTC  Collisions: 


Zij=^(av)ljdt=^(*v)ljdt 

Pij  =  w  (av)ij  dt  IV ceu 
Pmax  =  W  (av)™  dt/Vcell 

N select  =  ^Fniav)^  dt/Vcell 


Collide  if: 

pnnrl/1  \  ^  N collide  —  Pij 

Rand(1)  <  JU  -  P, ^ 


<CTv>r 


J 


a  - 


:  □ 


O  c*s  O' 


Review  of  Fractional  Collisions 


NTC  Collisions: 


o  (Collision  Rate  Volume): (Cell  Volume) 
o  Select  Fraction  of  \ N 2  Possible 
o  Probability  of  Event  Ratio  of  Volumes 
o  Correct  Non-Equilibrium  Frequency 

Fractional-NTC  Collisions: 


Zij=^(av)ljdt=^(*v)ljdt 

Pij  =  w  (av)ij  dt  IV ceu 
Pmax  =  W  (av)™  dt/Vcell 

N select  =  ^Fniav)™  dt/Vcell 


Collide  if: 

pnnrl/1  \  ^  N collide  —  Pij 

Rand(1)  <  JU  -  P, ^ 


<CTv>r 


a  - 


4 


:  □ 


O  c*s  O' 


Review  of  Fractional  Collisions 


NTC  Collisions: 


o  (Collision  Rate  Volume): (Cell  Volume) 
o  Select  Fraction  of  \ N 2  Possible 
o  Probability  of  Event  Ratio  of  Volumes 
o  Correct  Non-Equilibrium  Frequency 

Fractional-NTC  Collisions: 


Zij=^(av)ljdt=^(*v)ljdt 

Pij  =  w  (av)ij  dt  IV ceu 
Pmax  =  W  (av)™  dt/Vcell 

N select  =  ^Fniav)™  dt/Vcell 


Collide  if: 

P  n  n  HC  I  ^  N collide  — 

Rand(1)  <  JU  -  P, ^ 


<^v>r 


o  Select  /  by  Cost/Accuracy  Tradeoff 


A" select  —  f  Ap 


S  > 


4 


:  □ 


O  c*s  O' 


Review  of  Fractional  Collisions 


NTC  Collisions: 

o  (Collision  Rate  Volume): (Cell  Volume) 
o  Select  Fraction  of  \ N 2  Possible 
o  Probability  of  Event  Ratio  of  Volumes 
o  Correct  Non-Equilibrium  Frequency 

Fractional-NTC  Collisions: 


Zij=^(av)ljdt=^(*v)ljdt 

Pij  =  w  (av)ij  dt  IV ceu 
Pmax  =  W  (av)™  dt/Vcell 

N select  =  ^Fniav)™  dt/Vcell 


Collide  if: 

P  n  n  HC  I  ^  N collide  — 

Rand(1)  <  JU  -  P, ^ 


<^v>r 


o  Select  /  by  Cost/Accuracy  Tradeoff 
a  Collision  Aw  Scaled  for  Skipped 


A" select  —  f  Ap 

= S£z» 


S  > 


:  □ 


O  c*s  O' 


Review  of  Fractional  Collisions 


NTC  Collisions: 


o  (Collision  Rate  Volume): (Cell  Volume) 
o  Select  Fraction  of  \ N 2  Possible 
o  Probability  of  Event  Ratio  of  Volumes 
o  Correct  Non-Equilibrium  Frequency 

Fractional-NTC  Collisions: 


Zij=^(av)ljdt=^(*v)ljdt 

Pij  =  w  (av)ij  dt  IV ceu 
Pmax  =  W  (av)™  dt/Vcell 

N select  =  ^Fniav)™  dt/Vcell 


Collide  if: 

P  n  n  HC  I  ^  N collide  — 

Rand(1)  <  JU  -  P, ^ 


<av>r 


N select  f  Np 

=  S£z» 


Wi=Wi  -  A Wij  &  Wj=Wj  —  Awy 

W’(W„+l)=AWy  &  w{At,+2)=Aw,j 


&  > 


o  Select  /  by  Cost/Accuracy  Tradeoff 
o  Collision  Aw  Scaled  for  Skipped 
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SWPM  Collisions 


Stochastic  Weighted  Particle  Method: 


o  Developed  by  Rjasanow  &  Wagner 
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Review  of  Conservati 


Merge  to  Pair  — »  DOF  for  Conservation: 

o  (n+2):2  yields  Exact  Mass, 
Momentum,  and  Kinetic  Energy 
Conservation 

o  Applied  Spatially  also  Shown  to 
Conserve  Electrostatic  Energy 

o  Though  Energy  Conserving, 

Still  Thermalizes  VDF 


Merge 


_  ^(rc+2)  ^ 

Wi  I  Vi  - 


Wcell  =  Er  ’  W‘ 

‘  _  1  y^(«+2) 

WCe//  2— // 

1  Y^(W+2) 


V2  =  ^  E;"+2  V- (Vi  -  A 

WCe//  1  \  1  ) 

W(a/b )  =  Wm/2 
V(a/b)  =  V  ±  Tl\fV1 

Similarly:  x,(a/b)  =  *  ± 
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Collisional  Beams  in 


©  Initial  Bi-Max wellian  Distribution  in 
Potential  Well 

©  NTC  Collisions  Results  in  Beam 
Thermalization 

©  Fractional-NTC  Collisions  Produce 
Same  Behavior 

o  Particles/Cell  Dramatically  Different 

o  Fringe  Extends  to  Lower  Densities  with 
Variable  Weights 


itential  Well 
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-0.5  0  0.5 
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Number  Density,  logl 0(  n  ) 
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Collisional  Beams  in 


o  Initial  Bi-Max wellian  Distribution  in 
Potential  Well 

®  NTC  Collisions  Results  in  Beam 
Thermalization 

®  Fractional-NTC  Collisions  Produce 
Same  Behavior 

o  Particles/Cell  Dramatically  Different 

o  Fringe  Extends  to  Lower  Densities  with 
Variable  Weights 

o  Relative  ‘Error’  Unknown  without 
Analytical  Solution  or  High  Fidelity 
Simulation 
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Mach  2  Argon  Shock 


ID  Normal  Argon  Shock  Test 
o  Simple  Verification  vs.  DS IV 

o  Initial  Conditions: 

To  =  293K,  n0  =  lE22/m3,  v0  =  637.4(m/s) 
o  Initial  Jump  to  Post-Shock  at  1cm 

o  VHS  Collisions: 

7Ve/=273K,  dre/=4.17A,  ujvhs={ 0.81 
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Mach  2  Argon  Shock 


ID  Normal  Argon  Shock  Test 

®  Simple  Verification  vs.  DS IV 

o  Initial  Conditions: 

To  =  293K,  n0  =  lE22/m3,  v0  =  637.4(m/s) 
o  Initial  Jump  to  Post-Shock  at  1cm 

©  VHS  Collisions: 

Tref =21 3K,  dref=A.\l A.,  ccvz/.s— 0.8 1 

©  Time  Average: 
n  from  t£  [80, 100)/xs 
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Normalized  Density 


TURF  -  SWPM+Octree 


Target  N/Cell  Quadrupled  per  Line 
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Mach  2  Argon  Shock 


ID  Normal  Argon  Shock  Test 

®  Simple  Verification  vs.  DS IV 

o  Initial  Conditions: 

To  =  293K,  n0  =  lE22/m3,  v0  =  637.4(m/s) 
o  Initial  Jump  to  Post-Shock  at  1cm 

©  VHS  Collisions: 

Tref =21 3K,  dref=A.\l A.,  ccvz/.s— 0.8 1 

©  Time  Average: 
n  from  t£  [80, 100)/xs 

©  Error  (Normalized  Li): 

err=\n  —  n\/n 

©  Error  Controlled:  err  oc  y/N /cell 
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TURF  -  SWPM+Octree 


Target  N/Cell  Quadrupled  per  Line 
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Mach  8  Argon  Bow  Shock 


2D  Argon  Shock  Test 

©  Initial  Conditions  like  M=2  Except: 
vo  —  2550m/s 

©  Specular:  x=5  —  5.04mm  with  )-±2mm 

©  Half  Domain  Modeled: 

80/xm  x  80/xm  Cells 


TURF:  n  -  Standard  DSMC 


X  (mm) 

TURF:  n  -  SWPM+Octree 


Mach  8  Argon  Bow  Shock 


2D  Argon  Shock  Test 

©  Initial  Conditions  like  M=2  Except: 
vo  —  2550m/s 

©  Specular:  x=5  —  5.04mm  with  )-±2mm 

©  Half  Domain  Modeled: 

80/xm  x  80/xm  Cells 

©  Time  Average: 
n  from  tG  [80,  100)/xs 

o  SWPM  Similar  to  Standard  DSMC 


TURF:  n  -  Standard  DSMC 


X  (mm) 

TURF:  n  -  SWPM+Octree 
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Mach  8  Argon  Bow 


2D  Argon  Shock  Test 

©  Initial  Conditions  like  M=2  Except: 

vo  —  2550m/s 

©  Specular:  x=5  —  5.04mm  with  y=±2mm 

©  Half  Domain  Modeled: 

80/xm  x  80/xm  Cells 

©  Time  Average: 
n  from  te  [80,  100)/xs 

o  SWPM  Similar  to  Standard  DSMC 
o  Despite  Different  Np /  Cell 
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Distribution  A:  Approvei 
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TURF:  n  -  Standard  DSMC 


X  (mm) 

TURF:  n  -  SWPM+Octree 


TURF  Np/Cell  -  Standard  DSMC 


X  (mm) 


TURF  Np/Cell  -  SWPM+Octree 
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Issue  with  Collide  then  Merge 


o  Larger  Nseiect  Better  Approx,  of  Collision  Integral 
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Issue  with  Collide  then  Merge 


Larger  Nseiect  -»  Better  Approx,  of  Collision  Integral 
f-NTC  Produces  2x-Particles  per  Nseiect  mf  Np 
Particle  Memory  Requires  oc  Nmax  (1  +  2f)Nmax 


Issue  with  Collide  then  Merge 


o  Larger  Nseiect  Better  Approx,  of  Collision  Integral 
o  f-NTC  Produces  2x-Particles  per  Nseiect  =fNp 
o  Particle  Memory  Requires  oc  Nmax  (1  +  2f  )N,nax 
9  For  DSMC-like  Results,  /  «  0(1) 
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Issue  with  Collide  then  Merge 


o  Larger  Nseiect  -+  Better  Approx,  of  Collision  Integral 
o  f-NTC  Produces  2x-Particles  per  Nseiect  =fNp 
o  Particle  Memory  Requires  oc  Nmax  (1  +  2f  )N,nax 
9  For  DSMC-like  Results,  /  «  0(1) 

9  Time  Accurate  or  Dense  Simulations,/  «  0(10)+? 
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Issue  with  Collide  then  Merge 


9  Larger  Nseiect  -+  Better  Approx,  of  Collision  Integral 
9  f-NTC  Produces  2x-Particles  per  Nseiect  =fNp 
9  Particle  Memory  Requires  oc  Nmax  (1  +  2f  )N,nax 
9  For  DSMC-like  Results,  /  «  0(1) 

9  Time  Accurate  or  Dense  Simulations,/  «  0(10)+? 
9  Merge  Contracts  back  to  0(Nmax)  Particles 
9  Merge  Immediately  after  Collide  per  Spatial  Cell?.. 

9  Sort  for  Merge  still  oc  (1  +  2 f)  log (l  +  2/)? 


& 


□ 


'O  c*s  O' 


Issue  with  Collide  then  Merge 


9  Larger  Nseiect  -+  Better  Approx,  of  Collision  Integral 
9  f-NTC  Produces  2x-Particles  per  Nseiect  =fNp 
9  Particle  Memory  Requires  oc  Nmax  (1  +  2f  )N,nax 
9  For  DSMC-like  Results,  /  «  0(1) 

9  Time  Accurate  or  Dense  Simulations,/  «  0(10)+? 
9  Merge  Contracts  back  to  0(Nmax)  Particles 
9  Merge  Immediately  after  Collide  per  Spatial  Cell?.. 

9  Sort  for  Merge  still  oc  (1  +  2 f)  log (l  +  2/)? 

9  Combine  Collision  and  Merge  in  Single  Step? 
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Collide  to  Bins 


o  Fractional  Collision  as  Rate  Equation 
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Collide  to  Bins 


o  Fractional  Collision  as  Rate  Equation 
o  Bin  Moments  needed  for  Particle  Pairs 
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Collide  to  Bins 


o  Fractional  Collision  as  Rate  Equation 

o  Bin  Moments  needed  for  Particle  Pairs 

©  Particle  Pairs  (i,j)  Picked  Randomly 

o  DSMC-like  Collision  (VHS,VSS,etc.) 
Random  x,  0  -w  (v*/ ,  vy) 
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Collide  to  Bins 


o  Fractional  Collision  as  Rate  Equation 

o  Bin  Moments  needed  for  Particle  Pairs 

©  Particle  Pairs  (i,j)  Picked  Randomly 

o  DSMC-like  Collision  (VHS,VSS,etc.) 
Random  x,  0  (v*/ ,  vy) 

©  Octree  to  Find  i'  and  /  Bins 
8L  Few  Levels  to  Search 
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o  Fractional  Collision  as  Rate  Equation 

o  Bin  Moments  needed  for  Particle  Pairs 

©  Particle  Pairs  (i,j)  Picked  Randomly 

o  DSMC-like  Collision  (VHS,VSS,etc.) 
Random  x,  0  (v*/ ,  vy) 

©  Octree  to  Find  i'  and  /  Bins 
8L  Few  Levels  to  Search 

Conserve  Mass,  Momentum,  and  Energy 
Memory  Constant  Independent  of  Nsdecl 
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Mach  2  Argon 


ID  Normal  Argon  Shock  Test 
®  Mach  2  Case  Repeated 


Average(err) 


-  B2B 
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Mach  2  Argon  Shock  -  B2B 


ID  Normal  Argon  Shock  Test 
©  Mach  2  Case  Repeated 
©  Bin-to-Bin  Collsions  Results  Similar 


TURF  -  SWPM+Octree 


TURF  -  Bin  to  Bin 


Mach  2  Argon  Shock  - 


ID  Normal  Argon  Shock  Test 
©  Mach  2  Case  Repeated 
©  Bin-to-Bin  Collsions  Results  Similar 

©  Target  Np/  Cell  Still  Error  Control 
(Target  N/Cell  Quadrupled  per  Line) 
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X  (m) 


TURF  -  Bin  to  Bin 


Mach  2  Argon  Shock  - 


ID  Normal  Argon  Shock  Test 
©  Mach  2  Case  Repeated 
©  Bin-to-Bin  Collsions  Results  Similar 

©  Target  Np/  Cell  Still  Error  Control 
(Target  N/Cell  Quadrupled  per  Line) 

©  Collision  Core  ~3x  Slower 
©  Non-Ideal:  Dynamic  Range  Low 
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Average(err)  Average(err) 


TURF  -  Octree 


X  (m) 


TURF  -  Bin  to  Bin 


Mach  2  Argon  Shock  - 


ID  Normal  Argon  Shock  Test 
©  Mach  2  Case  Repeated 
©  Bin-to-Bin  Collsions  Results  Similar 

©  Target  Np/  Cell  Still  Error  Control 
(Target  N/Cell  Quadrupled  per  Line) 

©  Collision  Core  ~3x  Slower 
©  Non-Ideal:  Dynamic  Range  Low 
©  Proof-of-Concept  with  Real  X-Section 
o  Expansion/Plume  will  be  Better  Case 
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X  (m) 


TURF  -  Bin  to  Bin 


Mach  8  Argon  Bow  Shock 


2D  Argon  Shock  Test 

o  Mach  8  Case  Also  Repeated 


TURF:  n  -  Standard  DSMC 


X  (mm) 

TURF:  (n)  SWPM+Octree 
TURF:  Np/Cell  -  Standard  DSMC 


X  (mm) 
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Mach  8  Argon  Bow  Shock 


2D  Argon  Shock  Test 

o  Mach  8  Case  Also  Repeated 
o  Bin-to-Bin  Collsions  Results  Similar 
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Mach  8  Argon  Bow 


2D  Argon  Shock  Test 

©  Mach  8  Case  Also  Repeated 
©  Bin-to-Bin  Collsions  Results  Similar 
©  Target  Np/ Cell  Still  Error  Control 
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Distribution  A:  Approvei 


TURF:  Np/Cell  -  Standard  DSMC 


X  (mm) 


TURF:  (Np/Cell)  Bin-to-Bin 
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Mach  8  Argon  Bow 


2D  Argon  Shock  Test 

9  Mach  8  Case  Also  Repeated 
9  Bin-to-Bin  Collsions  Results  Similar 
9  Target  Np/ Cell  Still  Error  Control 

o  B2B  Run  with/=4x  Collisions 
(Note:  SWPM+Octree /=  lx) 


Standard  -  Collisions 

Standard  -  Total  Run 

548.9s 

7945.3s 

lx 

100% 

SWPM+Octree  -  Collisions 
SWPM+Octree  -  Total  Run 

2719.6s 

9542.4s 

4.95x 

120% 

Bin-to-Bin  -  Collisions 
Bin-to-Bin  -  Total  Run 

13163.6s 

18860.5s 

24.0x 

237% 

Robert  Martin  (AFRL/RQRS) 


Distribution  A:  Approvei 


TURF:  Np/Cell  -  Standard  DSMC 
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TURF:  (Np/Cell)  SWPM+Octree 
TURF:  Np/Cell  -  Standard  DSMC 


X  (mm) 


TURF:  (Np/Cell)  Bin-to-Bin 
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Mach  8  Argon  Bow 


2D  Argon  Shock  Test 

9  Mach  8  Case  Also  Repeated 
9  Bin-to-Bin  Collsions  Results  Similar 
9  Target  Np/ Cell  Still  Error  Control 

o  B2B  Run  with/=4x  Collisions 
(Note:  SWPM+Octree /=  lx) 


Standard  -  Collisions 

Standard  -  Total  Run 

548.9s 

7945.3s 

lx 

100% 

SWPM+Octree  -  Collisions 
SWPM+Octree  -  Total  Run 

2719.6s 

9542.4s 

4.95x 

120% 

Bin-to-Bin  -  Collisions 
Bin-to-Bin  -  Total  Run 

13163.6s 

18860.5s 

24.0x 

237% 

9  Some  Cost  Compensated  by  Lower  Np 
9  Too  much  Fill  for  Better  Wake 

9  Significant  Optimizations  Still  Needed 
(i.e.  Data  Structures,  Sort->Sums,  v-Bounds, 
Morton  curve) 
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TURF:  (Np/Cell)  SWPM+Octree 
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TURF:  (Np/Cell)  Bin-to-Bin 


vse;  Distribution  Unlimited;  PA  #16326 


16/18 


o  Standard  Collision  Incompatible  with  Variable  Weight 
©  SWPM+Octree  Option  for  Variable  Weight  Collision 
o  Bin-To-Bin  Potentially  Alleviates  Memory  Constraints 
o  Initial  Verification  vs.  Standard  Shock  Cases  Positive 
o  Limited  Utility  in  Standard  Shock  Cases 
©  Performance  with  Strong  Expansion/Plume  Needed 
©  SWPM/Bin-to-Bin  more  Useful  for  Trace  Species? 
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